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SOME EFFECTS OF FREQUENCY ON THE CONTRIBUTION OF А. VERTICAL ТАП, TO THE 
FREE AERODYNAMIC DAMPING OF A MODEL OSCILLATING IN YAW! 
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SUMMARY 


The directional damping and stability of a fuselage—vertical- 
tail model oscillating freely in yaw were measured at a Mach 
number of 0.14 and compared with the damping and stability 
obtained by consideration of the effects of unsteady lift. The 
contribution of the vertical tail to the damping % yaw of the 
model agreed unthin the limits of experimental accuracy unth 
predictions of the approximate finite-aspect-ratio unsteady-lift 
theories for the range of frequencies tested, and the two- 
dimensional unsteady-lift theory predicted a much greater loss 
in damping with reduction in frequency than was shown by 
experiment or by the approximate finite-aspect-ratio unsteady-lrft 
theories. For frequencies comparable to those of lateral airplane 
motions, both the unsteady-lift theories and the experimental 
results indicated that the contribution of the vertical tail to the 
directional stability of the model was relatively independent of 
frequency. 

INTRODUCTION 


The advent of high-speed airplanes of high relative density 
has focused attention on certain problems associated with 
the dynamic stability of airplanes which, because of previous 
unimportance, have heretofore been neglected. Among these 
problems is the effect of the periodicity of the airplane motion 
on the effective values of the various stability derivatives. 
Reference 1 and, more recently, references 2 and 3 have 
indicated the possibility of sizeable effects from this problem. 

An appreciable amount of theoretical work on these 
unsteady-lift effects exists at present; however, only a small 
amount of experimental substantiation of the results is 
available for the low-frequency range of oscillation.. As a 
result, а program has been undertaken in the Langley sta- 
bility tunnel to determine the effects of such variables as 
frequency and amplitude of motion on the contribution of 
the various airplane components to the stability derivatives 


of present-day airplane configurations. The work reported, 


herein covers that phase of the investigation which considers 
frequency effects on the directional damping and stability 
of а model undergoing a freely damped oscillatory yawing 
motion. The effects of vertical-tail aspect ratio and сош- 
pressibility as predicted by the theoretical treatments are 
discussed in relation to the experimental stability charac- 
teristics obtained by the free-oscillation and by the curved- 
flow procedures. 


SYMBOLS 


The data are referred to the system of stability axes and 
are presented in the form of standard NACA coefficients of 
forces and moments: about the quarter-chord point of the 
mean aerodynamic chord of the normal wing location of the 
model tested. (See fig. 1.) The coefficients and symbols 
used herein are defined as follows: 


C; lift coefficient, 1/45, 

Cp drag coefficient, 2/49. 

Ст lateral-force coefficient, Y /qS., 

С, rolling-moment coefficient, 1/|49.0 
Oa pitching-moment coefficient, M/qS ut w 
CO, 01 yawing-moment coefficient, /09,0. 
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Relative wind 


2222 


2 
Section А-А 
FigURB 1.—System of stability axes. Arrows indicate positive forces, 
moments, and angular displacements. Yaw reference generally 
chosen to coincide with initial relative wind. 


1 Supersedes МАОА TN 2657, “Bome Effects of Frequency on the Contribution ofa Vertical Tail to the Free Aerodynamic Damping of a Model Oscillating in Yaw” by John D. Bird, 


Lewis В. Fisher, and Sadie M. Hubbard, 1952 
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lift, Ib 

drag, Ib 

lateral force, Ib 

rolling moment, ft-lb 
pitching moment, ft-lb 
yawing moment, ft-lb 


dynamic pressure, = p УЗ, lb/sq ft 


mass density of air, slugs/cu ft 
free-stream velocity, ft/sec 
area, sq ft 

span, ft 

chord, ft 

mean tail chord, ft 


5/3 
mean aerodynamic chord, a f edy, ft 
. 0 


spanwise distance from the plane of 
symmetry, ft 

aspect ratio, 02/5 

angle of attack, deg 

angle of sideslip, radians 


| 


time, sec 


yawing-velocity parameter referred to 
semispan of wing 

angle of yaw, radians 

yawing velocity, d [dt, radians/sec 

y&wing acceleration, 020/012, radians/sec? 

angular velocity, radians/sec 

frequency, cps 

reduced-frequency parameter гакі to 
semichord of vertical tail, wc,/2V 

mechanical spring constant, ft-Ib/radian 

yawing moment of inertia, ft-Ib-sec? 

time to damp to one-half amplitude, sec 

damping constant . 

period of oscillation, Sec 

nondimensional tail length хен to, semi- 


chord of vertical tail, — 
с 3 2 


distance from origin of axes to midchord 
point of vertical tail, ft 


circulation functions used by Theodorsen 


(ref. 6) 

finite-span functions used by Biot and 
Boehnlein (ref. 9) 

finite-span correction used by Reissner and 
Stevens (ref. 8) 

compressible-flow functions used in place of 
Е and 6 (refs. 10 and 11) 


A,B functions of k defined in text 
_ об, 
Е. Tí | 
2V 
9C, 
QE: 
д ұя 
д 
Or 
РАЙ 
90, 
ng Op 
| оС, 
Cu Bb 
o Ре 
2V 
oC. 
Org og. 
` ON 
oN 
МУ;=— 
па Я 
ON. 
Мв 
oN 
N; SS 
к В 
Subscripts: 
w wing 
i vertical tail 
0 amplitude 
7 due to friction 
APPARATUS 


MODEL AND OSCILLATING STRUT 


The model used in this investigation and shown in figure 
2 had separable wing, fuselage, and tail surfaces and was 
constructed almost entirely of balsa wood in order to mini- 
mize the mass and make the natural frequency of the model 
on its mounting as high as possible. The configuration is 
representative of that of several present-day high-speed 
airplanes. The aspect ratio А, and mean chord с, of the 
vertical tail were 2 and 4.9 inches, respectively. 

The model was mounted at the quarter chord of the wing 
mean aerodynamic chord on a l-inch-diameter rod which 
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33.60 


Figure 2.—Drawing of model tested. АП dimensions are in inches. 


rotated with respect to a circular hollow strut and which 
was connected to it by means of flexure pivots made of 
Swedish steel. Flexure pivots were used rather than bear- 
ings in order to minimize the friction in the system. The 
pivots were alined with the center of the rod on which the 
model was mounted. The entire assembly, which is shown 
in figures 3 and 4, was fastened to the ceiling of the tunnel 
test section and covered by а streamlined fairing. | 

The inner rod or oscillating strut, which was permitted to 
rotate and upon which the model was mounted, extended 
through the ceiling of the tunnel test section and was capped 
by а horizontal crosspiece. The period of oscillation of the 
model was controlled by varying the moment of inertia 


= до ттоо ор 


Tunnel ceiEng 


Flexure pivots 
-Fairing я 
Fixed strut” с 
` Wind 
МӨР 
СЭНС? 
| „-Моде! 
МІГ”. tunnel È 


FraunE 3.—Diagram of test equipment. 


. about the vertical axis of the system. This variation was 


accomplished by using crosspieces of different lengths and 
by clamping weights at the ends of the crosspieces. 


RECORDING DEVICE 


A continuous record of the displacement in yaw of the 
model was made on photographic film by means of an optical 
recording system. A mirror clamped to the section of the 
oscillating strut which extended outside the tunnel reflected 
a beam of light from a point light source within the recorder. 
The beam was focused on a roll of photographic film by 
means of an 8-inch-focal-length double-concave lens. The 
film roll was driven at constant speed by gearing it to an 
electric motor. A 10-cycle-per-second neon timing light 
incorporated in the recorder exposed timing lines on the film 
in order that time as well as model displacement could be 
measured accurately on the film record. 


TESTS 
OSCILLATION TESTS 
The free-damping tests were conducted with the model in 
ап inverted position, in order to minimize strut interference 


effects, and at zero angle of attack with respect to the 
fuselage reference line. The tests were conducted on the 
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Кісовв 4.—Test model mounted on free-oscillation strut in Langley 
stability tunnel. 


combination of fuselage and tail surface &nd on the fuselage 
alone because preliminary tests had shown that the presence 
of the wing had no apparent effect on the damping at zero 
angle of attack. The range of moments of inertia and fre- 
quencies covered are shown together with the test results 
in table I. Results of runs made with wind on and off are 


given in order to illustrate the magnitude of the tares 


applied. 
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When tests were made of the fuselage alone, a torque rod 
was employed to provide a stable restoring moment to the 
fuselage which was equal to the aerodynamic moment pro- 
vided by the vertical tail when attached. The torque rod 
was clamped between the extended oscillating strut and a 
building structural member as shown in figure 3. The 
spring constant of the system was 6.8 foot-pounds per radian 
for the flexure pivots alone and 32.8 foot-pounds per radian 
for the combination of the flexure pivots and the torque rod. 


FORCE TESTS 


The steady-state stability characteristics of the model 
were determined by means of wind-tunnel force tests through 
an angle-of-attack range from 0° to 8° for the fuselage alone, 
the fuselage and tail surfaces, and the complete configura- 
tion. The model was mounted on a single-strut support and 
all forces and moments were measured by a six-component 
balance system. The steady-state sideslipping derivatives 
were determined through a sideslip range of +5°, and the 
yawing derivatives were determined at values of the yawing- 
velocity parameter 7b/2V of 0, —0.0291, —0.0616, and 


- —0.0810 by standard Langley stability-tunnel curved-flow 


testing procedures. 

All tests were conducted in the 6- by 6-foot test section of 
the ‘Langley stability tunnel at a dynamic pressure of 24.9 
pounds per square foot, а Mach number of about 0.14, and 
а Reynolds number of 442,000 based on the wing mean 
aerodynamic chord. Standard jet-boundary corrections 
have been applied to the angle-of-attack, the drag-coeflicient 
and the rolling-moment-coefficient data obtained from the 
force tests. In addition, the steady-state yawing data have 
been corrected for the effects of the lateral static-pressure 
gradient peculiar to the curved-flow testing procedure. 
Jet-boundary corrections were not applied to the oscillation 
tests because of their insignificance for the test conditions 
involved. See reference 4. 


REDUCTION OF OSCILLATION TEST DATA 


From the continuous film record taken of the motion of the 
model after an initial displacement, the amplitudes of the 


TABLE I 
EXPERIMENTAL RESULTS . 


Fuselage and tail 
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successive cycles were measured and plotted to a logarithmic 
scale against time. Inasmuch as the damping is of a loga- 
rithmic nature, the resulting plot is a straight line from which 
may be read the time for the motion to damp to one-half 
amplitude (бэ). 

For a system having a single degree of freedom in yaw 
such that for the stability-axis system the angle of sideslip 8 
is the negative of the angle of yaw y, the equation of motion 
can be expressed as 


U;—Nj3)$—(N,—N;4N,)0—(N,—N34v—0 (1) 


where Му, and Ny refer to the mechanical friction of the 


system and the flexure-pivot spring constant, respectively. 
The solution of equation (1) can be written as 


p=e-™(A sin 2nft+B сов 2xft) (2) 
where 
ММ [Му МЕХ», 
Iz—N} 2(Iz— N3) 
N,—N$--N;, 
2(Iz—N 3) 


m=— 


The damping constant may be expressed in terms of the time 
to damp to one-half amplitude as 


_ 0.693 
bisa 


The directional stability and damping of the system may 
then be expressed in nondimensional form in a manner similar 
to that of reference 5 as 


(bey 
Ont OCE) = 


Е ==МҮЗЎ | eee Е 
Cnr Өз; Б 49.0, 1» ал (61/9) tna orf йуз ч) 


The terms С», and C,; are necessarily in cothbination with 
С, й and C, , respectively, because of the nature of the motion 
being considered. For the frequencies employed in these 
tests, the factor C,.k? (22) should be small in comparison 
with Спи: 

In determining (й рз) „гна ог, Which corrects for the damping 
due to friction in the system, the model was oscillated at 
zero wind velocity with an equivalent mass in-place of the 
vertical tail and with a torque rod of sufficient strength to 
produce the proper frequency inserted in the system. The 


influence of the т Му 


2 др Us И (3) 


= term is very small compared with that of 


the inertia Мак, i omisión (4) and was neglected in these 
tests. 

The damping derivatives were obtained from data similar 
to the sample records of model motion and plots of logarithm 
of amplitude against time given as figure 5. А small residual 
yawing motion of the model attributable to turbulence in the 
tunnel airstream necessitated determining the time to damp 
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to one-half amplitude by the slope of the amplitude envelope 
at the largest angles of yaw in order to minimize the error in 
determining the damping derivatives. This source is 
believed to have resulted in no more than a 5-percent 
uncertainty in the results. 


ANALYTICAL CONSIDERATIONS 


. Various unsteady-lift theories were employed to calculate 
the contributions of the vertical tail to the directional 
stability and daniping of the model for comparison with the 
experimental results. ‘These theories included the incom- 
pressible cases covered by Theodorsen (ref. 6) for two 
dimensions and by Jones (ref. 7), Reissner and Stevens (ref. 
8), and Biot and Boehnlein (ref. 9) with various degrees of 
approximation for finite span. The two-dimensional work 
including the subsonic effects of compressibility as given by 
Possio and Frazer and Skan is also employed. A short 
summary indicating procedures for obtaining the necessary de- 
rivatives from these theories is included here for completeness. 


TWO-DIMENSIONAL DERIVATIVES 


The work of Theodorsen on the derivation of the expres- 
sions for the unsteady lift and moment of a two-dimensional 
surface undergoing sinusoidal oscillations in an incompres- 
sible fluid is given in reference 6. From this source, the 
moment per unit span about the airplane center of gravity 
of a vertical tail mounted on an airplane which is performing 
sideslipping oscillations can be shown to be 


МВ) = - rqe? 


8 att) свв] Ө 


A similar expression for ‘an airplane performing yawing 
oscillations is 
NW=—§ че 2—00 1)» ўе ita уы (6) 


An appendix to this report gives a discussion of: the con- 
version of these results from reference 6. 


(тава га mes 
ш ШІЛІГІ Мі ІТ АПА. ТПО " | | m ШЇЇ ЇЇ 


(a) Film record of yawing oscillation. 
(b) Logarithmic plots of amplitude variation with time. 
Етатвв 5.—Representative film record and resulting logarithmic plots. 
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The Theodorsen function C(k) is a complex circulation 
function of the reduced frequency and in terms of its real 
and imaginary parts is written | 


СФ F--iG 


The functions Е and G may be computed from expressions 
given in reference 6 in terms of Bessel functions of the first 
and second kinds of orders zero and one and of argument 
k. In stability practice equations (5) and (6) are usually 
interpreted in terms of derivatives as follows: 


N(8)— C, a 8,b.8--C,, 297 
S. e. (7) 
N()—C., 4 Зада а, Api V 
If the sinusoidal motions 
B= Ве“: 
y= yer 


are employed for а substitution into equations (5), (6), and 
(7), the real and imaginary terms may be equated and the 
derivatives Сер Сар С, ‚апа Ca; determined. The results of 


such an analysis would be as follows: 
кеде be v (a--5) Е 
5,0.) 1506) т] 
На 
==, s. (5; 2 |- EH e) 1-7)| 


The velocity terms contain only the F term of the Theodorsen 
function, whereas the acceleration terms contain the G term 
divided by Е. The corresponding side-force derivatives are 
derived in the appendix. For the free-oscillation tests 
conducted for this report, the motion is such that the angle 
of sideslip is the negative of the angle of yaw. The damping 
derivative measured therefore is the difference between the 
derivatives C,, апа С, з determined from equation (8). 
This fact was illustrated previously with respect to the 
reduction of the oscillation test data. In a like manner, the 
directional stability parameter measured can be shown to be 
a combination of the C, в and Ca; terms, specifically 


Combination of the results of equations (8) to obtain the 
factors calculated for comparison with the oscillation tests 
gives 


(8) 


25 В 
S,k 


Cn gt Са H X ig 8, 4 


С,,- Си аи 
(9) 
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where 


B=—(a~5 k—(a+5) 28 (à) ФЕ 


=(=+$ 5) (4+; 5) + (e-i —1) эра 


FINITE-SPAN DERIVATIVES 


The aerodynamic-span effect is considered in reference 7 
for wings of aspect ratios of 6 and 3 by correcting the aero- 
dynamic inertia and the angle of attack of the infinite-span 
surface. An approximation employed in this analysis 
concerning motions of long wave length may make the 
results subject to question for values of k near zero. The 
reference presents expressions for calculating the finite-span 
values of the Theodorsen function C(k). The use of these 
Е and С circulation functions in equations (8) and (0) 
results in values of the various derivatives considered for 
aspect ratios of 6 and 3. 

The aerodynamic-span éffect is considered for wings of 
arbitrary aspect ratio in reference 8. The three-dimensional 
effect of the finite span may be obtained by adding a correc- 
tion term с to the basic Theodorsen two-dimensional func- 
tion. The finite-span function C(k)+ø is employed in a 
rather lengthy and tedious computation in order to deter- 
mine the components of the unsteady lift for each of five 
points along the span of the oscillating surface. The various 
stability derivatives may then be determined by a graphical 
integration of the in-phase and out-of-phase components. 
An approximation made in this analysis is of such a nature 
as to make the method primarily suited to the higher aspect 
ratios. A practical lower limit of applicability is not known. 

А one-point approximation to Ше span effect for the cal- 
culation of the unsteady lift may be employed rather than 
the five-point procedure. Comparative calculations made 
for the present investigation have shown the one-point 
approximation to yield results for low aspect ratios which 
are essentially the same as those calculated by the longer 
procedure. 

Biot and Boehnlein (ref. 9) considered the aerodynamic- 
span effect on unsteady lift by an approach that yields a 
closed form for a one-point approximation to the span effect. 
The finite-span circulation functions used are 


P=F+iG 
Q—H-HJ 
For infinite aspect ratio 
| P=Q=C(k) 


where C(k) is the two-dimensional Theodorsen function. 


The functions P and Q may be calculated for the argument 
k by methods given in reference 9, and equations (8) and 
(9) are again used to calculate the various stability deriva- 
tives. The factors 4 and В which appear in equations (9) 


must, however, be written in terms of P and 0 as follows: 


A=(2a—NakG+(a—5) ы/а) 124-20F-+H 
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В--(2а-- NakF-+(a—5) kH—(a—5) k—2aG—J 
EFFECT OF COMPRESSIBILITY 


Subsonic compressibility effects on the unsteady damping 
in yaw and directional stability of a two-dimensional oscil- 
lating lifting surface can be considered through the use of 
reference 10 which is an extension of the calculations of 
Possio (ref. 11). The infinite-span compressible-flow values 
for the stability derivatives can be obtained by the use of 
equations (8) and (9) for Mach numbers 0, 0.5, and 0.7. In 
the notation of Frazer and Skan (ref. 10), however, 


«н (м+2-32)+м, 


B=4M,— 5(5+а) (м--2-227-м. 

The correspondence between the notation of Theodorsen 
and that of Frazer and Skan is shown in reference 12. The 
values of the 2 and М functions are tabulated by Frazer and 
Skan for Mach numbers 0, 0.5, and 0.7 and a significant 
range of reduced frequencies. Reissner has recently ex- 
tended his work on oscillating wings of finite span to include, 
to the same degree of approximation as for the incompressible 
саво, the subsonic effects of compressibility. (See ref. 13.) 


RESULTS AND DISCUSSION 


PRESENTATION OF RESULTS 


The experimental values of the longitudinal, directional, 
and yawing characteristics of the complete model, the model 
with wing removed, and the fuselage alone as determined by 
force tests are given in figure 6. The experimental values of 
the damping in yaw of the fuselage and tail assembly 
(designated FT) and the fuselage alone (designated Е) 
obtained by the freely damped oscillation technique for 
several frequencies of oscillation are shown in figures 7 and 8. 
The frequency and aspect-ratio effects on the damping in 
yaw contributed by the vertical tail as calculated from the 
unsteady-lift theories discussed previously are also shown in 
these figures. Figure 9 gives a comparison of the relative 
importance of C,, and C,; on the demping in yaw of the 
model. Figures 10 to 12 are similar figures summarizing 
frequency and aspect-ratio effects on directional stability. 
Figures 13 and 14 give the computed effects of compressi- 
bility on the vertical-tail contribution to the damping in yaw 
and directional stability of the test model. 


DISCUSSION 


The variation of the experimentally determined damping 
in yaw of the test model with frequency is in reasonably 
good agreement for the range covered by the tests with the 
various approximate theoretical treatments which consider 
the effect of finite span. (See figs. 7 and 8.) The contribu- 
tion of the fuselage to the damping in yaw is small for the 
range of test frequencies. A two-dimensional approxima- 
tion to the damping as is frequently employed for flutter 
work appears to be remarkably good for k>0.1 and for aspect 
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ratios of 6 and greater (Яр. 8). А large effect of aspect ratio, 
however, is shown for aspect ratios below 6 for the entire 
frequency range and for all aspect ratios for frequencies less 
than k=0.1. The two-dimensional result indicates a loss in 
damping at low frequencies much greater than the finite- 
span result and, in fact, shows a reversal in sign of the 
damping at sufficiently low frequencies. A number of papers 
dealing with unsteady lift have noted this characteristic of 
the two-dimensional aerodynamic forces. (See refs.. 1 to 3 
and 14.) 

An examination of the values of the derivatives О, and 
Оһ, ав computed from reference 6 for infinite aspect ‘ratio 


id from reference 7 for aspect ratios of 3 and 6 (fig. 9) 
indicates that the derivative Саҙ is critically dependent on 
reduced frequency for the low es of reduced frequencies 
and infinite aspect ratio. The derivative C, is only slightly 
dependent on the reduced frequency for infinite aspect ratio 
and even less so for the lower aspect ratios. The dependence 
of Съ, on reduced frequency is responsible for the decrease 
in damping shown for the infinite-aspect-ratio case as the 
frequency is reduced (figs. 7 and 8). 

The aspect-ratio-2 results shown in figures 7 and 10 were 


obtained from the theory of Jones (ref. 7) for aspect ratio 


3 by use of & correction factor which was the ratio of the 
lifting-line-theory lift-curve slopes for aspect ratios of 2 and 
3. This correction is only approximate and was made in 
order to place the results on a comparable basis. The varia- 
tion of C, — Оһ) with reduced frequency, shown in figure 7, 
indicates no appreciable effect of frequency to zero frequency. 
In fact, the result of Jones for aspect ratio 3 reduces to a 
finite value at k=0 given by 


C, — Са (й st) те (1.247--0.370а--0.515) 


which is not the case for the two-dimensional result that 
approaches infinity as Е approaches zero. This result is not 
expected to be very accurate close to £—0, however, because 
of the nature of the assumptions made in its derivation as 
mentioned previously. 

When the finite-span effects are ЖЕНТ) it appears that 
these effects are especially significant with regard to the 
damping in yaw of this model for aspect ratios smaller than 
6. (See fig. 8.) The theory also indicates a significant 
decrease in the damping at low frequencies when the aspect 
ratio was increased from 20 to o. 

The experimental values of damping in yaw obtained by 
free oscillation agree closely with the values of С, obtained 
by the standard curved-flow testing procedure used in the 
Langley stability tunnel. The curved-flow technique results 
іп the k=0 value of О,. No unsteady effects, however, 
are obtained by this technique. 

The unsteady-lift theories indicate that, at frequencies of 
oscillation from k=0 to k=0.1, the frequency has little 
effect on the directional stability contributed by the vertical 
tail but that high frequencies of oscillation can cause the 


directional stability contributed by the vertical tail to de- 


crease greatly and eventually to change sign (figs. 10 and 11). 
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Angle of attock, a, deg 
Fiaurs 6.—The steady-state stability characteristics of the model tested. 
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ЕТІ Curved-flow tests 
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FigunzE 7.—The effect of frequency on the damping in yaw for a low range of the reduced frequency. Comparison of several methods for cal- 
' culating frequenoy effect on the damping іп yaw. 
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Етсовв 8.—The effect of aspect ratio on the барва in yaw as calculated d the method of Biot and Boehnlein for a range of reduced frequencles. 
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(b) Range of reduced frequencies. 
FicurE 9.—The effect of reduced frequency on the stability derivatives 
Са, and Caz 


Changes in aspect ratio, for low aspect ratios, have a large 
affect on the directional stability; for high aspect ratios, 
shanges in aspect ratio are less significant. The experimental 
values of directional stability obtained from the free-damping 
ests remained constant for the range of frequencies investi- 
rated, although the magnitudes are somewhat smaller than 
hose predicted for the tail contribution by the calculation 
nethods for aspect ratio 2. This result is not unexpected 
зесарве all the finite-span unsteady-lift theories employed 
naintain some resemblance to the Prandtl lifting-line 
heory which, of course, yields values of the lift-curve slope 
эг low aspect ratios which are unrealistically high. The 
alues of directional stability obtained by both the free- 
scillation and the static-test procedures were about the 
ime magnitude. 

Figure 12 gives a comparison of the contribution of a 
ertical tail to the stability derivatives Cs, and Cy, as 
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computed from references 6 and 7. These derivatives were 


combined in the proportion 


mu 
снощи (се) 


for comparison with the results of the oscillation tests. 
Although С. is of appreciable magnitude itself, the factor 


Ам ; 
С.Е? (2%) is very small relative to C,, for the range of 


frequencies around Е--0.01 and thus should not greatly 
affect the directional stability of the usual airplane. 

The computed two-dimensional Mach number effects on 
the damping in yaw and directional stability of the tail 
surface of the test model are shown in figures 13 and 14. 
These results were included to illustrate the powerful effect 
of compressibility on the damping contributed by the tail 
surface. A complete coverage of one-degree-of-freedom 
flutter by use of these aerodynamic forces which may be 
obtained from references 10 to 12 is given in reference 14. 
Compressibility effects decrease the tail-surfate damping at 
low frequencies and increase the damping at higher frequen- 
cies. "The directional stability increases at all frequencies 
with increase in Mach number with the largest increases 
occurring in the high-frequency range of oscillations. Refer- 
ence 15 which presents the two-dimensional oscillating 
aerodynamic forces at sonic speed indicates larger deviations 
from the results for incompressible flow than for the subsonic 
results given here. 


CONCLUSIONS 


The following conclusions were drawn from the results of 
an investigation made at a Mach number of 0.14 of the 


‘effects of frequency and aspect ratio on the directional 


damping and stability of a fuselage and vertical tail oscil- 
lating freely in yaw: 

1. The contribution of the vertical tail to the damping in 
yaw of the model agreed within the limits of experimental 
accuracy with predictions of the approximate finite-aspect- 
ratio unsteady-lift theories for the range of frequencies 
tested. 

2. Two-dimensional unsteady-lift theory predicted a much 
greater loss in damping with reduction in frequency than was 
shown by experiment or by the approximate finite-aspect- 
ratio unsteady-lift theories. 

3. Both the unsteady-lift theories and the experimental 
results indicated that the contribution of the vertical tail to 
the directional stability of the model was relatively independ- 
ent of frequency for frequencies comparable to those of 
lateral airplane motions. 


LANGLEY ÁERONAUTICAL LABORATORY, 


NATIONAL ADVISORY COMMITTEE FOR ÁERONAUTICS, 
Lawanzy Ево, Va., November 26, 1951. 
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Ғтаонв 11.—The effect of aspect ratio on the directional stability as calculated by the method of Biot and Boehnlein 
for a range of reduced frequencies. 
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FrcunE 12.—The effect of reduced frequency on the stability deriva- 
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(а) Lower range of reduced frequencies. 
(b) Range of reduced frequencies. 
Егаове 14.—The effect of compressibility on the directional stability. 


APPENDIX 


CONVERSION OF UNSTEADY FORCES AND MOMENTS FROM FLUTTER TO STABILITY AXES 


The conversion of the contribution of the vertical tail to 
the yawing moment of an airplane in terms of unsteady-lift 
functions and of the Theodorsen flutter axes to the stability 
axes which are fixed in the body is included here as a matter 
of general interest. Reference 16 gives considerable atten- 
tion to such axis conversions. The following sketches indi- 
cate the paths and attitudes of airplanes performing sinusoidal 
sideslipping and yawing oscillations within the limits of the 
stability-axis system: 


Sideslipping, B= 861% 


Yawing, уро" 
The Theodorsen coordinates А and а and the stability-axis 


coordinates В and ў are shown with respect to the paths 
followed by the airplanes in sideslipping and yawing flight 
In sideslipping flight 


— - eee —— — ——— m — ма ММ” 


— á— MX — ——— - 
= 
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and in yawing flight 
= -5 h=a | 


From reference 6 the contribution of the vertical tail to the 
yawing moment of an airplane in terms of the flutter coordi- 
nates is Е 


Ма, t= — 28 « (5-% Ус, re (1 gt +a?) a 502 i + 


34 


If this moment is for а sideslipping motion or a yawing 
motion in the sense given by stability axes fixed in the body 
and if the following terms are substituted for the sideslipping 
case 


pV с T (2+5) С(Е) | Vatht% 


I 


and for the yawing case 


= 
== 
the following results are obtained for the sideslipping case 


рс? (ажс У N, p Veer 
м) (FE д) е (at 
and for the yawing сазе 


N@=—2 i G- DATE (Gta) Ў Ў |+ 


пе (н) oG) i 


which reduce for the sideslipping case to 


5) СС 


М(В)=—тас жор зу Нача) сев | 
and for the yawing case to 


н - 2 


Nw=—§ ае? Ў 5-009 G- 


Similar expressions may be obtained for the side force in pure 


sideshpping and yawing flight. Reference 6 gives the side 
force of the vertical tail in terms of the flutter axes as 


У (о,ћ)= ы (Упарта. 
трус [Гаа 3-4) 3 


which, when the following terms are substituted for the 
sideslipping case 


and for the yawing case 


аттұ/ 
а=} 


acad 


becomes for pure sideslipping and yawing flight 
үф) а: | с, peu 


Y ()esgqos 2 z| 200 (5- ЗЕ ta % 


Substituting the motions, іп the manner done in the text, 
for the yawing moment, 


and 


B= poe?! 
p= poor 


into these expressions and into the expressions in terms ot 
stability derivatives which are 


| S00 1 
Y (8)=Or,8 gS wt Or,6 124" 


2 
Ү (у)= Оу} Bates Oyj Bere A 
gives the following relationships for the tail contribution tc 


the side-force derivatives when the real and imaginary part. 
of these equations are equated: 


S, C: 1 
Оу, ==2-=- S, bp TÉ 7) | 


Os - S (Ft) [0-2 2-4) 
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